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^ ■ Abstract 

We provide Standard Model expectations for the rare radiative decays B — > K*j, 
B — > p7 and B — > uj^, and the electroweak penguin decays B — ► K*£ + £~ and 
-B — > pl + l~ at the next-to-leading order (NLO), extending our previous results 
to 6 — ► (i transitions. We consider branching fractions, isospin asymmetries and 
direct CP asymmetries. For the electroweak penguin decays, the lepton-invariant 
mass spectrum and forward-backward asymmetry is also included. Radiative and 
electroweak penguin transitions in b — > d are mainly interesting in the search for 
new flavour-changing neutral current interactions, but in addition the B — > p7 
decays provide constraints on the CKM parameters (p,fj). The potential impact 
of these constraints is discussed. 



1 Introduction 



The radiative and electroweak penguin transitions b — > D7 and 6 — > D£ + £~ (D = d, s) 
provide valuable insight into the nature of flavour- changing neutral currents. Induced 
through quantum fluctuations in the Standard Model, they may be strongly influenced 
by new heavy particles. This has been studied extensively for inclusive b — > 57 decays 
over the past decade. Some years ago the QCD factorization approach to exclusive 
B decays [JJ E] was extended to exclusive radiative [SI IH Ej and electroweak penguin 
decays [3], thus opening the possibility to perform detailed studies of isospin breaking 
and CP asymmetries in these decays, as well as to obtain more precise predictions for 
branching fractions and forward-backward asymmetries. 

At the time when these papers were written, only the exclusive B — > K*^ decays 
had been observed, and a next-to-leading order (NLO) analysis could be done only for 
B — > i^*7, p7 and B — > K*£ + £~, the case of B — > p£ + £~ being excluded by the absence 
of the NLO (2-loop) virtual correction to the b — *> d£ + £~ transition. This missing piece 
of input has recently been computed [SJ Cj, while the formal justification of the QCD 
factorization method in radiative decays has been pursued in jHJ |U] • Moreover, the B 
factories have now put limits on the B — > p7 and B — > wf branching fractions [TU1 ITT] , 
and have performed first measurements of B — > K*£ + £~ [T21 ITU] including the lepton 
invariant mass spectrum and forward-backward asymmetry [Tlj . These theoretical and 
experimental advances motivate the following study of the observables of interest with 
a complete next-to- leading order calculation (except for weak annihilation). 

The paper is organized as follows: In Section El we give a brief introduction to the 
theoretical formalism. We specify the input parameters to the analysis and give the 
various contributions to the decay amplitudes in numerical form. The theoretical expec- 
tations for the relevant observables are summarized in Section El For both, B — > V'-f 
and B — > V£ + £~ decays (V = K*,p,u), we discuss branching fractions (lepton invariant 
mass spectrum and partially integrated branching fractions for B — > V£ + £~), isospin 
asymmetries (difference between charged and neutral B meson decay) and direct CP 
asymmetries. For B — > V£ + £~ the forward-backward asymmetry is also included. In 
Section 13.51 we collect the constraints on the CKM unitarity triangle that can be ob- 
tained from the branching fractions, isospin and CP asymmetries in B — > pj decays. 
The analysis of the radiative decays overlaps with recent work of Ali et al. [15 a and 
Bosch and Buchalla ^E] > and we compare our results to theirs in the appropriate places. 
We conclude in Section 0] The technical Appendix El contains the new decay amplitudes 
for B — > p £ + £~ related to the up-quark sector of the effective weak Hamiltonian. 

2 Theoretical input 
2.1 Formalism 

The formalism is described in some detail for b — > s transitions in jH] , and the extension 
to b — > d is straightforward. In the Standard Model the effective Hamiltonian for b — » d 
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Table 1: Wilson coefficients at the scale fi = 4.6 GeV in leading-logarithmic 
(LL) and next-to-leading-logarithmic order (NLL). Input parameters are 
= 0.220 GeV, m t (m t ) = 170 GeV, M w = 80.4 GeV and sin 2 ^ = 0.23. 
3-loop running is used for a s . 





Ci 


c 2 


c 3 


c 4 


c 5 


c 6 


LL 
NLL 


-0.5135 
-0.3026 


1.0260 
1.0081 


-0.0051 
-0.0048 


-0.0693 
-0.0836 


0.0005 
0.0003 


0.0010 
0.0009 




Cf 




Cg 


Cio 




°10 


LL 
NLL 


-0.3150 
-0.3094 


-0.1495 
-0.1695 


2.0072 
4.1802 




-4.3810 


4.2978 


-4.4300 



transitions can be written as 
with \W = V* d V qb , and 

10 

H® = c 1 o c l + c 2 o c 2 + ^c l o lJ 

2=3 

H$ = C x [0\ - Of) + C 2 [0% - O u 2 ). (2) 

For b — > s transitions is replaced by A^ „, so that the term X^H^ is CKM-suppressed 
and can be neglected in practice. The main technical point of the present paper is to 
add the decay amplitudes from H^g , which are relevant to the b — » d case. We have 
written the effective Hamiltonian in a form such that involves the differences of 
"tree" operators (cb) (Dc) and (ub) {Du) for b — > D (D = d, s) transitions. We use the 
operator basis as given in [3] . The numerical values for the Wilson coefficients at fi = rrib 
at leading-logarithmic (LL) and NLL order are collected in Table ^ The next-to-next- 
to-leading logarithmic (NNLL) results for Cg t \o were obtained from the solution to the 
renormalization group equations given in 0, including the recent computations of the 
3-loop mixing of the four-quark operators into Cg [17] and among themselves ^H]- The 
coefficient Cg is now complete at NNLL as formally required by our analysis. 

In the QCD factorization formalism the hadronic matrix elements are computed in 
terms of B meson form factors and hadron light-cone distribution amplitudes at leading 
power in a 1/mb expansion. Since the semi-leptonic operators C^io are bilinear in the 
quark fields, their matrix elements can be expressed directly through B — » V form 
factors. The other operators contribute to the decay amplitude only through the coupling 
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^(A« + Ai^)+ h - c - 



(1) 
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to a virtual photon, which then decays into the lepton pair. We therefore introduce 



(i*{q,riV(i/,e')\H$\B(p)) 



4tt 2 



-2iTl\q 2 )[EM B e^ - (e* ■ q)j/») 



-iT®(q*)(e*-q) 



Ml 



(3) 



where \V) denotes \p ) (\K* )) for B meson decay, and -V2\p°) or y/2\u°) (\K* )) for 
B° decay. In the heavy quark limit this matrix element depends on only two independent 
functions Tjf 1 corresponding to a transversely (a =_L) and longitudinally (a =||) polarized 
V. Most of the functions T® can be directly inferred from the calculation for the case 
B — > K*£ + £~ jH] with obvious replacements for K* — ► p(u). In Appendix El we give the 
new functions and point out the differences in for the p (&) meson compared 
to K*. The decay amplitudes can be further expressed as 
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p + — 
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(4) 



The second term incorporates hard scattering of the spectator quark. f B and 3>b,± refer 
to the S meson decay constant and light-cone distribution amplitudes, f\\ = f, f± and 
<p a to the corresponding quantities for light mesons. Furthermore = 1, H11 = my/E. 
The first "form factor" term is expressed in terms of the transverse and longitudinal 
"soft" form factors £j_ and £y, and C a , T a> ± denote perturbative hard scattering kernels. 
The significance of these terms will be discussed subsequently, but we mention here that 
in this work we use a definition of £y that differs slightly from [HUES] as explained below. 
In the description of B — > V7, there is no advantage in using the "soft" form factors. In 



this case we write T^' = Ti(0) C±'' + . . .. The coefficients C±' and T| v ^_ are related to 
the unprimed coefficients by (jHJ) given below. 



i(0' 



(i)> 



2.2 Input parameters 

A detailed discussion of the input parameters can be found in [3j. A summary is given in 
Table El where we have taken the CKM parameters from [20J - Unless stated otherwise, 
mj denotes the potential-subtracted (PS) heavy quark mass |2I] (see also the appendix). 
For the top quark mass we use the MS definition. 

The numerical value of Xb,+, related to the first inverse moment of the B meson 
light-cone distribution, is taken from the QCD sum rule calculation [22] • The decay 
constants and Gegenbauer moments of the light meson distribution amplitudes follow 
the values given in Table 1 of (231- We assume large uncertainties for these parameters, 
which cover in particular the recent evaluations of the first Gegenbauer moment of the 
kaon [21JI2H] 1 (renormalization-scale dependent quantities are evaluated at scale 2 GeV). 

-"^In our convention the momentum fraction used in the light-cone distribution amplitudes refers to 
the outgoing quark. Hence our a\(K*) is opposite in sign compared to these papers. 
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Table 2: Summary of input parameters and estimated uncertainties. 



M w 80.4 GeV 
m t (m t ) 170 ± 5 GeV 

«em 1/137 

Aq£v 5) 220 ± 40 MeV 
m 6iPS (2GeV) 4.6 ±0.1 GeV 
m c 1.5 ±0.2 GeV 
f B 200 ± 30 MeV 


AB,+(1.5GeV) 0.485 ±0.115 GeV 
t b o, t b ± 1.54 ps, 1.67 ps 
A 0.2265ig;oo23 

a omitS 

P 0.1891SZ 
f? 0.358+S-SS 


f K * >± 175 ± 25 MeV 
/*-.,[! 218 ± 4 MeV 
Oi(^*)j.,|| 0.2 ±0.2 
a 2 {K*)^ )]l 0.1 ±0.3 


/ pM ,i 150 ± 25 MeV 

ip,|| , / w ,|| 209 ± 1 MeV, 187 ± 3 MeV 

a>i(p,uj)_ L) \\ 

a 2 (>o)x, ||, ^H^H 0.1 ±0.3, 0.0 ±0.3 


M B tK;\\{0)/(2m K *) 0.47 ±0.09 
^^(O) 1 0.36 ±0.07 


M B Z Pj \\(0)/(2m p ) 0.37 ±0.06 
C P) ±(0) 0.27 ±0.05 



^In Section mi we determine £r:*,-l(0) = 0.26 from experimental data. This value rather than the 
one in the Table is then used in the subsequent analysis. 



We find it convenient to slightly change the convention for the longitudinal "soft" form 
factor compared to ^] . Denoting by a tilde the form factors in the old convention, 
the "soft" form factors are defined by the relations 

M B + m v 

«ll(0 2 ) = A,tf) Uf) = M<?) - A 2{q % (5) 

where A||(g 2 ), given by Eq. (66) of [3], determines the perturbative correction to the 
form factor relations between A , A\ and A 2 . Since A||(0) = 1 and A||(g 2 ) = 1 at tree 
level, the only difference is a change in the g 2 -dependence of £||(<? 2 ) versus i\\{q 2 )- 

The g 2 -dependence of the QCD form factors V, A\ and A 2 has been parameterized 
in the form 

F(n 2 ) = m 

{q 1 {l-a F q*/Ml + b F q±/M%y {) 

We use the definitions © and this equation with parameters given in Table |3] to deter- 
mine the g 2 -dependence of £j_ and £n. The normalization of the longitudinal form factor 
at q 2 = is 

{| <0) = ^<0) = ^*<0)-^*<0), (7) 
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Table 3: Form factor parameters from QCD sum rules [211] ■ The value V(0) 
is calculated from 11(0) via (jHJ) and (jSJ). 



meson 


V(0) 


ay 


by 


A 1 (0) 


ai 




A 2 (0) 


a 2 


b 2 


T\(0) 


P 


0.313 


1.37 


0.315 


0.261 


0.29 


-0.415 


0.223 


0.93 


-0.092 


0.285 


K* 


0.424 


1.55 


0.575 


0.337 


0.60 


-0.023 


0.283 


1.18 


0.281 


0.379 



but we do not use V(0) to obtain £j_(0). We use instead the relation between V(0) and 
Ti(0) (at scale /x = mj) [1^ to write 

ri (0)-^(0)^l-— j+— ^ U (8) 

with /x/ ~ 1.5 GeV, and obtain £j_(0) from the value of T"i(0) given in This ensures 
that, since the radiative decays B — > involve only the tensor form factor Ti or, 
alternatively, we obtain identical results for these decays independent of whether 
we express the decay amplitude in terms of the QCD form factor T\ or the "soft" form 
factor. 

Except for the longitudinal decay constant and second Gegenbauer moment, the 
hadronic parameters of the cu-meson are assumed to be identical to those of the p-meson 
for lack of more information. As a consequence we do not have meaningful estimates of 
the difference between the p° and u decay rates. 



2.3 Decay rates and distributions 

The decay amplitudes for B — > V'y and B — > V£ + £~ can be written in terms of [3 

Mi) _ Tj. \0) _ cit n eR , 



q 2 U<7 2 ) ' 



where i — t,u refers to the two different CKM factors. 2 



2 There is a factor of Am (q 2 ) missing in Eq. (41) of [3]. With the new definition of the £|| this factor 
is absent from the definition of . Furthermore, we now define by dividing through the QCD 

tensor form factor (at scale [i = mb) rather than £±. Hence C7 is expressed in terms of the coefficients 
Cf' and T { 1 ] ' + defined at the end of Section IOC see also the appendix). 
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With these definitions, we obtain the decay rate for B — > p{ui)^ in the form 

r(^^ h ) = gMls(i-gJ^^r l( o,^ 

x {lAl"'^"! 2 + \\^C^\ 2 - 2|A<'"Ai'"|Re (e"C^C? *) } , (10) 

with 5 = 1 for p~, and S = 1/2 for p° and u;. my denotes the mass of the light meson, 
which we include in the phase space factor, while in general m v terms are neglected. The 
dominant term is jA^C^I 2 ~ | | 2 , but the interference term is non-negligible and can 
be the source of interesting CP-violating and isospin-breaking effects. The CP-conjugate 
B decay follows from (|1U|) by the replacement a — > —a. 

The decay B — > p(uj)£ + £~ has a richer kinematic structure. Defining q 2 , the invariant 
mass of the lepton pair, and 6, the angle between the positively charged lepton and the 
B meson in the center-of-mass frame of the lepton pair, and summing over final state 
polarizations, the decay information is contained in the double differential distributions 




The lepton mass is set to zero, so this result applies to £ = e,/i. The CP-conjugate B 
decay follows from by the replacement a — > —a. The terms with angular dependence 
(l±cos 2 6) correspond to the decay into transversely and longitudinally polarized p(o;)'s, 
respectively. The term proportional to cos# generates a forward-backward asymmetry 
with respect to the plane perpendicular to the lepton momentum in the center-of-mass 
frame of the lepton pair. 
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Table 4: Breakdown of the decay amplitudes C 7 for B — > p7 and ory. 





/»(*) 
c 7 




Cf 
C(i)' 


-0.309 
-0.060 - 0.013a 




0.063 + 0.055a 


T(°) (annih.) 

T (i)/ (p) 
( W ) 




-0.009 - 0.012a 
-0.009 - 0.013a 




-0.031 -0.012* 
-0.029 -0.013* 


a® /nib (annih.) 


U.UUy \p ) 
-0.005 (p°) 
-0.004 (w) 


n IOC / r%~ \ 

— U.izo [p ) 
-0.012 (p°) 
0.011 (u) 


(P-) 

«!M (p°) 
(«) 


0.001 + 0.000a 
0.000 - 0.001a 
0.001 - 0.001a 


0.000 - 0.004a' 
0.000 + 0.002a- 
0.000 + 0.002a- 


(p") 

sum (p°) 

H 


-0.369 - 0.024a' 
-0.383 - 0.026a' 
-0.383 - 0.027* 


-0.093 + 0.039* 
0.019 + 0.045*- 
0.045 + 0.042a' 



2.4 Overview of amplitudes 

We briefly discuss the qualitative features of the decay amplitudes from which the main 
characteristics of the observables that we analyze in the following section can be deduced. 
Beginning with B — > V'-f we see from Tabled that the process is dominated by the short- 
distance electromagnetic penguin amplitude proportional to C| . However, there is an 
important radiative correction C^' to the quark matrix element ( , yd\H^ u ' \b) J2Z], as well 
as a sizeable correction from spectator scattering [HJ E] specific to the exclusive decays. 
It is worth noting that these two corrections add constructively in the top-sector of the 
effective Hamiltonian (C^p), but tend to cancel in the up-sector. The most distinctive 
feature of exclusive h —> transitions is the large weak annihilation contribution to 
B^ — ► p ± 7, which has been discussed extensively jU I2H I2H EO]- Although formally 
suppressed by a factor of l/m&, the annihilation amplitude is the largest contribution 
to df for B ± — > p ± 7, because it is enhanced by a Wilson coefficient ten times larger 
than Cy . The leading annihilation amplitude is short- distance dominated and can be 
computed in the heavy quark limit. However, as discussed in the Appendix, there is a 
large theoretical uncertainty associated with this computation due to both parameter 
uncertainties (/#, Xb+) and power corrections. The remaining power-suppressed ampli- 
tudes listed in Tableware small effects and relevant only to the extent that they provide 
the leading source of difference between p° and u, or, in the case of b — * s transitions 
K* + and K*° j3*T| . The Table shows that the difference of p° and uj is negligible, but 
it should be remembered that the hadronic parameters of the u have been set equal to 



7 



those of p° in the absence of better information. 

From a phenomenological point of view, the complete information about the four 
different B — > pj decays is encoded in the parameters |(C7*' ) ) p0 | ~ 0.38, 

{C ' ))P+ 1 = 5 -±±3± ~ -0.04 + 0.00*, (13) 



and 



(M u hp° (C^) p+ 

1 7 J = e e ie ° ~ -0.06-0. 1H, 1 7 J =e + e ie + ~ 0.24 - 0.12i. (14) 
(C?)p° (C?V + 

The discussion of the theoretical errors of these parameters is deferred to the following 
sections. The decay amplitudes are then proportional to 1 — R ut eie % ^ i±a ^ (i = 0,+), 
where 

X^/xf = -R ut e ia (15) 

with R ut — 0.46, and a ~ 94° in the Standard Model, and the upper (lower) sign refers 
to the decay of a B (B) meson. 

The amplitude structure of B — > V£ + £~ is somewhat more complicated due to the 
presence of the axial- vector short-distance contribution proportional to C w , the duplica- 
tion of amplitudes for transverse and longitudinal polarization of the vector meson, and 
the g 2 -dependence of these amplitudes. For very small q 2 the amplitude is dominated 
by the photon pole and exhibits a behaviour qualitatively very similar to that discussed 
above for a real photon. In most of the region, in which our theoretical treatment is 
applicable (q 2 < 7GeV 2 ), C w and the longitudinal amplitude determine the decay char- 
acteristics with the exception of the forward-backward asymmetry. The latter is directly 
proportional to the expression in square brackets in the last line of (jlljl . and is therefore 



sensitive to the real part of Cg^ and e ±,a Cg"|, which is expected to go through zero in 
the g 2 -range of interest. Table El gives the values for Cg„ (a =_L, ||) at q 2 = 5 GeV 2 for 
the decays to the p meson. We note again that the two-loop correction to the quark 
matrix element {^f*d\H^ \b) is rather large, in fact as large as the contributions from 
four-quark operators at leading order in a s contained in the functions Y^^q 2 ). This is 
true in particular for the amplitude in the up-sector, which was calculated only recently 
jH] (see also [?]). Similarly to the radiative decay the power-suppressed weak annihilation 
amplitude is the largest remaining term in the transverse amplitude of the decay to p ± . 
Turning to the longitudinal amplitudes C^'^\ we should recall that we do not include 
l/m^-corrections in this case. The reason for this is that the important isospin-breaking 
effects exist already at leading power as is evident from the entry for (annih.) in 
Table EJ This comes from a weak annihilation contribution to B — > p£ + £~, which is 
not suppressed by factors of a s or l/m^ relative to the short- distance amplitude and 
which now also has a non- vanishing (and rather uncertain) absorptive part. This leads 
to a large difference of the longitudinal amplitude in the up-sector between the charged 
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Table 5: Breakdown of the decay amplitudes Cg^ (q 2 = 5 GeV 2 ) for B — > p£ + £ . We 
define aCf = (2m b M B /q 2 ) Cf for C$? ± and aCf = [2m b /M B ) Cf for C^j. 





At) 


u 9,± 


r (t) 

9,|| 


p(u) 
9.11 


Cg 

F(g 2 ) 
aCf 

cw 


4.30 
0.51 + 0.06i 
-3.01 
-0.82 - 0.14i 




-0.07- 0.84i 


0.23 + 1.27z 


4.30 
0.51 + 0.06i 
-0.54 
-0.39-0.01i 




-0.07- 0.84i 


-0.07 + 0.72i 


(annih.) (p°) 
T« (p-) 




-0.21 - 0.15i 




-0.09-0.15i 


0.04- 0.03i 
-0.02 + 0.01i 
0.05 - 0.04i 
0.02 - 0.05i 


1.31-0.86i 
0.13-0.09? 
-0.08 - 0.07i 
-0.09 - 0.04i 




0.07- 0.03i 
-0.04 + 0.02i 


0.39 - 1.06i 
0.04-0. Hi 








0.03 + 0.00i 
-0.02 + 0.00i 


0.00 - 0.03i 
0.00 + 0.02i 






snm {p0) 


0.87- 0.25i 
0.72 - 0.2H 


0.46-0.81z 
0.10 + 0.20i 


3.97- 0.00i 
3.88 + 0.02i 


1.08- 1.05z 
-0.ll-0.25z 



and neutral p final state similar to the situation for Q (now including the absorptive 
part). The amplitude in the up- sector is correspondingly uncertain (see the discussion 
in the appendix). 

In the following sections we perform a detailed analysis of the observables of interest, 
including theoretical error estimates not represented in Tables 0] and El The tables should 
provide useful information to understand the numerics and parameter dependences of 
the various observables. 



3 Phenomenological analysis 

In this section we discuss several observables related to rare radiative B decays which 
are currently studied at B factories, or will be measured at future high-luminosity exper- 
iments. The B — > Vj decays have been recently studied within the QCD factorization 
framework by Ali et al. |15j . and by Bosch and Buchalla |lfij . and in most aspects our 
analysis leads to similar results. In the case of B —>■ K*£ + £ decays we perform an update 
of our previous work [3]. Finally, we extend the discussion to B — > p{u)£ + £~ decays. 
Here the inclusion of the recently calculated two-loop correction to the b —>■ d^j* vertex 
from light-quark loops jUl U\ turns out to have a significant numerical impact on the 
various decay asymmetries. In the following discussion we focus on the sensitivity of 
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Table 6: Experimental results for the CP-averaged B — > V'y branching ratios 
(in units of 10~ 6 ). 





BABAR [10, 32 


BELLE [11, 33J 


CLEO [34J 


Average 


Br(B° -> K*° 7 ) 
Br(5+ -> iT* + 7) 


39.2 ±2.0 ±2.4 
38.7 ±2.8 ±2.6 


40.1 ±2.1 ±1.7 
42.5 ±3.1 ±2.4 


45.5l£j ± 3.4 
37.6±|:§ ± 2.8 


40.1 ±2.0 
40.3 ±2.6 


Br( J B° -> wy) 
Bi(B° -> p° 7 ) 
Br(E+ -> p+7) 


< 1.0 

< 0.4 

< 1.8 


< 0.8 

< 0.8 

< 2.2 


< 9.2 

< 17 

< 13 


< 0.8 

< 0.4 

< 1.8 



the radiative decays to the theoretical input parameters, including the CKM elements, 
short-distance Wilson coefficients, and hadronic parameters. 

3.1 B — > V^f branching fractions 

We first consider the B — > K*j decays. Using the central values and theoretical uncer- 
tainties for the input parameters in Table |2j we obtain for the branching fractions 

Br(2?° - = (^f) 2 (7.4l°j|^i^| had ) • ID" 5 , (16) 

Br(5+ - JT+ 7 ) = (^pj (7.4i°j|^| h J . 10~ 5 . (17) 

The main theoretical uncertainty comes from the B — > K* tensor form factor, the other 
hadronic and CKM parameter uncertainties being small as indicated. 3 Comparison with 
TableElshows that with the form factor from light-cone sum rules, the theoretical value is 
too high compared to data. Since the short-distance physics explored in inclusive b — > S7 
decays is well in line with the standard model expectation, we must conclude that either 
there are large power corrections not included in the computation, or the tensor form 
factor is smaller than 0.38. The existence of large power corrections beyond the known 
annihilation terms would invalidate the basic assumption of the factorization approach 
that the heavy quark expansion provides a sensible approximation. In the absence of any 
direct determinations of the form factor, we consider the reduction of the form factor a 
viable option, and use the branching fraction data to obtain |3] 

^(0)l cxp = 0.28 ± 0.02. (18) 

3 The theoretical error is computed from the parameter uncertainties given in Table [21 including 
the renormalization scale dependence, but no attempt is made to quantify the error from higher-order 
perturbative corrections and power corrections. This should always be kept in mind in the interpretation 
of theoretical errors. 
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This number will be taken as a reference value in the following. 

The uncertainty in the form factors is also a major obstruction to a clean interpreta- 
tion of the B — > p7 decays. The situation for K* may suggest that the B —>■ p tensor form 
factor should be reduced in equal proportion, in which case we obtain for the branching 
fractions averaged over B decay and the corresponding CP-conjugate B decay 



Br{B° p° 7 ) 



td I 



8.25-10 



-3 



0.21 



(5.0t£f) 



10" 



(19) 



Br(B- 



P + l) 



\V t 



td 



8.25-10 



-3 



0.21 



10.3 



-1.5N 

-1.2) 



lo- 



rn 



On the other hand, if the origin of the large B —>■ K* form factor in QCD sum rules is 
a misestimate of SU(3) breaking effects, Tf (0) ~ 0.29 is not excluded and the central 
values of the branching fractions are rescaled by almost a factor of 2. The branching 
fraction for B° — > wy equals Br(B° — > p 0/ ~f) within the given accuracy, and under the 
assumption that the form factors are the same. Comparison with the experimental limits 
on exclusive b — > d'y transitions in Table H favours the scenario with a small form factor 
when the standard range of \V t d\ is assumed. In the following, however, we continue to 
use the default value Tf (0) = 0.29 ± 0.04. 

It is often stated that the ratios of B — > p 7 to B —>■ K*^ branching fractions are 
better suited for a determination of \Vu\- This is based on the assumption that the ratio 
F = Tf (0)/T 1 K *(0) is better known than the form factors themselves. Unfortunately, 
given that we do not know for certain whether the current estimates of the form factors 
are affected by a normalization or SU(3) breaking problem, we must assume 0.75 < F < 
1.05 at least. 4 The ratios of CP-averaged branching fractions are determined by 



Br(B c 



A) 



Br(B° -> iT* 7) 

= V 

2 



If* 

2 



V, 



td 



V, 



Is 



[ 1 - 2R ut e cos a cos 9 + R 2 ut e^] 



V, 



td 



Vt 



Is 



jl — 2R ut cos a [- 



•'.on- !!;!!!; • <• • 



[0.0218:8?]} 



(21) 



and 



Bt(B+ -> p+7) 



Vt 



Id 



Vt 



Is 



{1 



2R, llf e_i_ cos a cos 9 A 



K ut e + 



V 



td 



V 



is 



{l-2R ut cos a [0.2418;$] + i& 



[0.0718;^]} 



(22) 



where a small phase space correction {M\ — m 2 p ) 3 /(Mg — m 2 Kt ) 3 ~ 1.02, and a SU(3) 
breaking correction from IC^/IC^^* — 1 = ±0.02 can be safely neglected. 5 The global 

4 The lower limit corresponds to the assumption that QCD sum rules predict the ratio of form factors 
correctly. The upper limit assumes that T-f (0) is predicted correctly and (0) is determined by data 
which yields l|18|) . 

5 With Tf(0) = 0.21 instead of 0.29, the central values of the numerical entries in (|21(l change from 
(-0.06,0.02) to (-0.01,0.01) and in (jHJ from (0.24,0.07) to (0.40,0.18). 
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CKM fit j2D] returns R u t = 0.46±0.06 and cosa = — 0.07±0.20, so the interference term 
is expected to be suppressed. Since the value of e cos 9 is rather small due to a partial 
cancellation of the perturbative corrections C^- 1 ' and T*- 1 -* (see Table HJ), the B° — > p°7 
(and to a slightly lesser extent also B° — ► decay is rather insensitive to i? n t and 
a. The neutral ratio thus provides a direct constraint on |Vt<f/Ks| for a given value of 
the ratio of tensor form factors F = Tf(0)/T 1 K *(0). For the charged ratio e + cos# + is 
dominated by the weak annihilation contribution, to which we assign a 50% error (see the 
discussion in the appendix). This error is by far the largest uncertainty in B + — > p + 7, 
and makes the charged ratio somewhat less useful to constrain |Vtd/Vt s |. 

At this point it is appropriate to compare our results eocos#o = —0.06 ± 0.06 and 
e + cos# + = 0.24 ± 0.18 to fSJ GUI- Bosch and Buchalla employ a formalism identical 
to ours and give 0.0 ±0.1 and 0.4 ± 0.4 for these parameters. We find agreement with 
their central values, when we adjust the parameter Xb,+ to the smaller value 350 MeV 
used by them. Ali et al. do not give eocos^o and e+cos#+ explicitly, but the curly 
brackets in (|21I22|) for which we obtain 1.00 ± 0.01 and 1.03 ± 0.06 correspond to their 
1 + AR= 1.09 ± 0.07 and 1.12 ± 0.10. 



3.2 B — > V£ + £ branching fractions 

The differential branching fractions for B — > K*£ + £~ decays have been discussed in detail 
in our previous analysis [Sj. The present work includes changes of input parameters 
(see Table EI), the now available NNLO anomalous dimensions relevant for the Wilson 
coefficient Cg [T7J as well as a modification of the treatment of form factors. In 
particular, the transverse form factor is now taken from data via (|18jl. and we no longer 
use the simple 1/E n (n = 2,3) model for the g 2 -dependence of the form factors. The 
combined effect of these changes is to decrease the branching fraction by about 30% 
relative to j3J, almost exclusively due to the change in the form factor input. 

The g 2 -spectrum rises sharply for small q 2 , where it is dominated by the photon pole 
and possibly "contaminated" by hadronic resonances. On the other hand the factor- 
ization approach is valid only for q 2 < 7 GeV 2 , away from the charm threshold. We 
therefore advocate that measurements of exclusive branching fractions be compared to 
the integral 

60 V Br(B+ - {^£§P) 2 (3-3^) ■ ID"'. (23) 

1 GeV 2 

of the CP-averaged spectrum. The corresponding integral for the neutral B meson decay 
is about 10% smaller. The largest hadronic uncertainty arises from the longitudinal form 
factor, which has therefore been scaled out, such that only the residual Ao-dependence 
is included in the error. 

Recently the B factories have presented their first results for (partially) integrated 
B — ► K*£ + £~ decay rates O H3] • Of particular interest to us is the second bin of the 
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Belle measurement ^3], which translates into 



8GeV 2 



r , n dBi(B -> K*£+£-) , , , , 7 

7 d( ? J2 = ( 4 ' 8 -i^lstat. ± 0.3| syst . ± 0.3| model ) • 10" 7 , (24) 



4GeV 2 



dq 



which we interpret as an average of the charged and neutral decay. The spectrum is nearly 
flat in this range, so we may compare one half of this number to our theoretical result 
integrated from q 2 = 4GeV 2 to 6GeV 2 . We obtain (1.2 ± 0.4) • 1(T 7 (now including 
the uncertainty in A$ (0)), which is only about one half of the central value of the 
experimental result. It will be interesting to see whether this difference is due to a 
detection bias, or rather points to a problem with the theoretical input. 

A new result of the present work is the NNLL (next-to-next-to-leading logarithmic) 
prediction for the B — ► p£ + £~ decay rate below the charm threshold. Integrating the 
spectrum of the CP-averaged decay as above, we obtain 

**y> = UgL ^Y!> y (4 ,« :S) . 10 -, (25) 

1 GeV 2 




y rfg 2 V 8.25-10- 3 0.50 ' " Ll 

lGeV 2 X 

In Figure ^ we show the spectrum rfBr(5° — ► p°£ + £~)/dq 2 . The qualitative features 
are similar to 5 — > K*£ + £~, namely the NLO correction is very important at small q 2 , 
where the transverse amplitude is large, but at large q 2 the spectrum is flat, and the 
NLO corrections from different sources cancel to produce a negligible net effect. 

By the time the B — > p£ + £~ decays can be measured, the CKM elements, as well 
as the short-distance Wilson coefficients Cg and Cio for b —* s£ + £~ transitions, will 
be known with high accuracy. The B — > p£ + £~ decays can provide additional insight 
into the structure of flavour- changing neutral current interactions in the context of new- 
physics scenarios with non-minimal flavour violation, where Cg and Cf Q can differ both, 
in magnitude and in phase, from their counterparts in the b — » s sector. To achieve the 
necessary theoretical accuracy, it is crucial to control the hadronic uncertainties. We may 
expect the form factors at small q 2 to be known with much better accuracy than today, 
presumably from lattice simulations. In this situation the factorization formalism should 
allow us to obtain stringent constraints on Cg and Cf Q from the neutral decay mode. 
On the other hand, the theoretical control of charged decays remains more challenging 
because of the leading-power annihilation effects whose q 2 behaviour strongly depends 
on the model for the 5-meson wave function. 6 



6 These effects are suppressed in the branching fraction in the Standard Model, because cos a is small. 
This suppression need not hold in extensions of the Standard Model such as discussed here. 
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Figure 1: CP-averaged differential branching ratio for B° — ► p°£ + £~ at NLO as a function 
of q 2 (solid line, in units of 10 _9 GeV~ 2 ). The light (yellow) band shows the total 
theoretical uncertainty. In the dark (green) band, the uncertainties related to the CKM 
parameters and the form factor Aq(0) are excluded. The dashed line shows the LO result. 

3.3 Isospin asymmetries 

The isospin asymmetry in B — > K*^ has been discussed in [^TJ. Our calculation gives 
(all rates averaged over the CP-conjugate decay) 



A(K^) - nB^K^)-T(B^K^) _ _028_ 

J) ~ T(B° - K^) + r(£?+ -> ^*+ 7 ) ~ If* (0) ^- 2 - 9j/ °' ^ 



which is consistent with the experimental number (3.6 ± 3.8)% deduced from Table El 
As already discussed, isospin breaking in radiative decays is a power-suppressed effect. 
Although the calculation is believed to capture the dominant effect, its theoretical status 
is less certain than the calculation of branching fractions and CP asymmetries. Our result 
is smaller than the result of the original calculation [3T] , because we evaluate the Wilson 
coefficients at a scale of order rrif, rather than (^AgcD) 1 ^ 2 - The motivation for this 
choice is that the four-quark operators factorize below the scale /i ~ m^, and the gluon 
exchange between all quark lines responsible for the renormalization group running of 
the Wilson coefficients is no longer relevant at smaller scales. Our result is larger than 
the one given in ^H], partially because in this paper the isospin-breaking hard-scattering 
corrections (denoted by "o^/mb" in Table HJ are not included, but primarily because the 
larger QCD sum rule value of the K* tensor form factor is used there. 

The extension of the isospin analysis to B — > K*£ + £~ has been discussed in . We 
refrain from giving an updated discussion here, but only mention that with the different 
choice of scale described above, we expect an even smaller isospin asymmetry than the 
estimate in [33], which implies a larger sensitivity to isospin- violating new physics effects. 

While the isospin asymmetry in B — > K*^ decays mainly probes the magnitude of 
penguin Wilson coefficients in weak annihilation, the isospin asymmetry for fq decays 
is sensitive to CKM parameters through the interference with a large tree annihilation 
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Figure 2: Isospin asymmetry A(p 7 ) as a function of the CKM angle a. The band 
displays the total theoretical uncertainty which is mainly due to weak annihilation. The 
vertical dashed lines limit the range of a obtained from the CKM unitarity triangle fit. 



amplitude with a different CP phase. We obtain for the isospin asymmetry (again an 
average over the CP-conjugate decay is understood) 



A(p7) 



P + l) 



2T(B° -> p° 7 ) 



A fi+ 5 - 4 l + 5 - 8 l W 
-^•O -4.2lcKM-5.6lhadJ /0 - 



(28) 



The largest hadronic error (4%) comes from the weak annihilation contribution to which 
we assign a 50% error. The uncertainty labeled "CKM" illustrates the sensitivity to the 
CKM parameters. Using the abbreviations (jl3H4|) . and neglecting terms quadratic in 
8 + , eo — 0.1 and rj, we can approximate the exact result for the isospin asymmetry by 



A(py) ~ 5 + — 2 R ut cos a (e + cos 9 + — e cos 6* ) + R 2 ut et t 



-0 07fi +0 - 042 

u - u ' °-0.049 



2R ut cos a [0.30±gig] + R 2 ut [0.07^] 



(29) 



with a relative error of less than 10%. 7 For a near 90° the isospin asymmetry is predicted 
to be small, and the numerical result is dominated by the isospin asymmetry in the 
top quark sector, 8 + . Our result (|28J) therefore differs from where the term S + is 
neglected. Eq. (|29~|l also shows that there is a large sensitivity to cos a, although with a 
large uncertainty related to the weak annihilation amplitude. This dependence is shown 
in Figure 121 

The isospin asymmetry in the differential decay spectrum for B — > p£ + £~ is defined 
in analogy to (|28|). In the limit q 2 — > its numerical value approaches ([29)1 . For invariant 
lepton-pair masses well above the photon pole at q 2 = the decay spectrum is dominated 
by the longitudinal rate, and therefore the isospin asymmetry mainly comes from the 
amplitudes Cyl with i = t,u. It can be seen from Table^Jthat the dominant effect comes 

7 With T[(0) = 0.21 instead of 0.29, the central values of the numerical entries in (|29|) change from 
(-0.076,0.30,0.07) to (-0.104,0.42,0.18). 
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Figure 3: Isospin asymmetry A(p£ + £~) as a function of q 2 . The solid (long-dashed) line 
shows the next-to-leading (leading) order result for a = 94°. The band represents the 
hadronic uncertainty. The two dashed lines give A(p£ + £~) for a = 24° (lower curve) 
and a = 164° (upper curve). 

from the leading-power weak annihilation contribution [3| to the charged decay mode 
B + — > p + £ + £r which, as already explained, has a large uncertainty. When cos a is not 
near zero, the asymmetry is generated mainly by Cg*) and the sign of the asymmetry is 
determined by — cos a just as for B — > pj. On the other hand, for small cos a (expected 
in the Standard Model), the asymmetry is small, comes mainly from the top quark 
sector, and its sign is opposite to In any case, the effect is diluted compared to 

B — > p7 because of the isospin-symmetric contribution of the Wilson coefficient Cio to 
the longitudinal amplitude. This situation is illustrated in Figure El where we show the 
predicted asymmetry for a = 24°, 164° and 94° (default). For the last value of a the 
theoretical uncertainty is also shown. The maximum of the isospin asymmetry around 
q 2 ~ 1 GeV 2 can be explained by the observation that the longitudinal amplitude for the 
charged decay becomes dominated by the leading-power weak annihilation contribution 
for q 2 — > 0, since As j _(g 2 ) -1 increases logarithmically. Hence the isospin asymmetry 
increases (coming from larger q 2 ) until the transverse decay amplitude wins over and the 
isospin asymmetry turns to the negative value f)29j) for q 2 = 0. A reliable prediction is 
currently possible only in the larger q 2 region shown in the Figure. 

As mentioned above the B — > p£ + £~ decays are mainly of interest in the context of 
scenarios where physics beyond the standard model modifies the electroweak penguin 
coefficients Cg, Cio in b — > d transitions. If 9 denotes the extra CP-phase of Cg, the 
isospin asymmetry depends, roughly speaking, on a c s = a— 6. Since a will be known, the 
sensitivity of the isospin asymmetry to a shown in Figure H3 translates into a sensitivity 
to 9. The theoretical uncertainty of the Standard Model reference value due to weak 
annihilation implies that only a significant new phase could be unambiguously detected. 
In this case a significant modification of the absolute value of Cg and Cio is also probable 
and might be directly observable in the lepton invariant mass spectra. 
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Figure 4: Direct CP asymmetry in B — > p°7 (solid), B + — > p + 7 (long-dashed) and 
5 — >■ u;7 (short-dashed) decay as a function of the CKM angle a. The band shows the 
theoretical uncertainty for the case of B — > p°7. Note that we display minus the CP 
asymmetry. 



3.4 Direct CP asymmetries 

The direct CP asymmetries for B —>■ p7 decays are given by 



o r[E ^A]-r[ff -.p°7] = 2i^eosinflosina 
( l>l '" ' r[S° - pQ 7 ] + r[S0 - p0 7 ] l-2i? ut e cos^ cos« + J R2 te 2' l " ' 



and an analogous equation for the charged B decay with (eo, 9q) — > (e+, #+). The values 
of e, cos 9i and have already been given in (|21I22|) . Here we also need 

e sinfl = -0.111^1, e + sin# + = -0.12±°;g|. (31) 

The largest theoretical uncertainties are from the residual renormalization scale depen- 
dence (±0.03) and the charm quark mass (±0.02). This reflects the fact that a next-to- 
leading order calculation of the branching fractions results in leading-order predictions 
of direct CP asymmetries, which are therefore more sensitive to unknown higher-order 
and power corrections. It is worth noting that the product e + sin 9 + is less dependent 
on the weak annihilation contribution than the individual factors e + and sin#+, since 
the annihilation amplitude has no strong phase relative to the leading electromagnetic 
penguin amplitude. It follows from (|31|) that, within uncertainties, the CP asymmetries 
in neutral and charged B — > p7 decays are of similar size 

^cp(A) = (-10.4^| CKM ^| had )%, 

A CP (p + l) = (-10.7^| CKM ^| had ) %. (32) 

The dominant dependence of the CP asymmetries on the CKM parameters is through 
Rut sin a. The corresponding constraint in the (p,fj) plane is discussed in Section I3~5l In 
FigureEJwe show the dependence of the direct CP asymmetries on the CKM angle a. The 
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Figure 5: Direct CP asymmetries in B° — > p°£ + £r (lower set of curves) and B + — > p + £ + £~ 
as a function of q 2 . The solid (dashed) curves show the next-to-leading order (leading 
order) result. The widths of the bands represent the hadronic uncertainty. Note that we 
display minus the CP asymmetry. 

asymmetries for B — > p°7 and B — > uj'j are indistinguishable within uncertainties. Our 
result is in agreement with [1J ^3 within theoretical uncertainties, though [15.] displays 
a slightly larger difference between the neutral and charged CP asymmetries. 

The direct CP asymmetry arises in B — > p £ + £~ decays from the interference between 
the C^l{q 2 ) and Cg a (q 2 ) amplitudes with different strong phases. From Table El we 
deduce that the largest contributions to strong phases come from the one-loop function 
Y(q 2 ), from the coefficient function and, in the case of the charged decay, from 
annihilation topologies which involve the B — > 7* form factors ^jT (q 2 ) (see appendix). 

The two-loop virtual correction calculated in [HI E] plays an essential role here, since 
it cancels a large part of the imaginary part of Y^ u \q 2 ) in Cg"[ and Cg) . For the charged 
decay B^ — > p ± t r £~ the residual relative strong phase is then dominated by the rather 
uncertain annihilation effects. For the neutral decay B° — > p°£ + £~, on the contrary, the 
strong phases are expected to be rather small, leading to a small CP asymmetry within 
the Standard Model. Our numerical prediction for the two decay modes is shown in 
Figure El The increase of — Aqp(p + £ + £~) near q 2 ~ 1 GeV 2 occurs for the same reason 
as discussed for the isospin asymmetry and is correspondingly uncertain. In extensions 
of the Standard Model a is again replaced by a e ff, but given the theoretical errors and 
given that the CP asymmetry is expected to be nearly maximal in the Standard Model, 
it will be difficult to disentangle a new phase unless the CP asymmetry in the charged 
decay mode is strongly suppressed. 




-A CP (p°£+£-) 

q 2 [Ge\ 2 ] 
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3.5 CKM constraints 



As seen in previous sections, measurements of the branching ratios, isospin- and CP- 
asymmetries in B — > pry decays are sensitive to the CKM elements and thus provide 
interesting constraints on the parameters p and fj, which define the apex of the unitarity 
triangle. In this section we summarize these constraints. 

At present, there exists only an upper experimental limit for the B —>■ fry branching 
fractions. Interestingly, our theoretical result (|T9"j) already saturates the experimental 
bound, and therefore detection of this decay is expected in the very near future. The 
existing upper limit on B — > p°7 translates into a useful bound on \ Vtd\, which is comple- 
mentary to the bound from the non-observation of B s B s -mixing. The key point is that 
the curly bracket in (|21|) is constrained to be close to 1. Adding the dominant parameter 
dependencies linearly and doubling the error on the input values of (p, fj) we find 

0.94 < 1 — 2R ut eo cos a cos e Q + Rl t e 2 < 1.05. (33) 

Using the data from Table |H1 this translates into the bound 



td 



ts 



1 ± 0.03 / 2Br(5° 



^Br(£ c 




(34) 



where we used the conservative range F > 0.7, which gives the weakest bound. A similar 
bound has been obtained in jTB]- 8 This limit already cuts into the range |V^/V t ,s| = 
0.204±8;8|| (CL=0.05) obtained from the standard fit 20 . This is displayed in Figure El 
where the area to the left of the solid black curve is excluded by (|34j). and the standard 
values of (p, fj) are shown as a point together with their errors as in Table [21 9 

No data currently exists for the isospin and direct CP asymmetries. To illustrate the 
possible constraints we assume that they are observed with a value that corresponds to 



our theoretical expectations (|29l32j) without experimental error. The theoretical uncer- 
tainty from the input parameters (excluding CKM parameters) then translates into the 
dashed (isospin asymmetry) and dash-dotted (direct CP asymmetry in B — > p°7) bands 
in the Figure. The constraint from the CP asymmetry in the charged decay 5 ± — > p ±r y 
is very similar to the neutral one in the vicinity of the standard (p, fj) range and not 
shown. The shape of these constraints can be understood from the relations 

R ut sin a = ^ R ut cos a = 1 - - — 1 / (35) 

(l-p) 2 + r] 2 (l-p) 2 + r] 2 

They imply to good approximation that the direct CP asymmetry requires (p, fj) to lie 
on a circle of radius l/(2a) with center (1, l/(2a)) where 

^cp(p°7) 



2eo sin 6 



(36) 



8 In |15j an average over p°,p^,oj is performed, which reduces the significance of the constraint due 
to the larger theoretical uncertainty in B — > p ± 7 and the weaker experimental upper limit on the 
and u: final states. 

9 The solid line becomes slightly inaccurate far away from the standard range, since the dependence 
of Ip^jl on (p, fj) is neglected. 
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Figure 6: Constraints on the unitarity triangle from the ratio of the CP-averaged B — > 
p°7 and B° —>■ K* 0/ y branching fractions (solid black), the isospin asymmetry in B — > p7 
(dashed) and the direct CP asymmetry in B° — > p°7 (dash-dotted). The area to the 
left of the black line is excluded by the experimental limit on B° — ► p°7. See text for 
explanations. 

which always passes through (1, 0). Similarly, the isospin asymmetry constrains (p, fj) to 
lie on a circle whose center is always at fj — and whose radius is determined from (J29|) . 
We see from Figure |H1 that the two asymmetry constraints intersect nearly orthogonally, 
the intersection region being shaded in grey. This area is further constrained by the 
limit on the B — > p7 branching fraction. The three observables together (and a similar 
constraint from the direct asymmetry in 5 ± — > p ± '~f) demonstrate that the B — > p7 
decays alone provide valuable independent information on CKM parameters. Further- 
more, if inconsistencies with the standard CKM fit appeared, this would point towards 
anomalous effects in b — > d'-f transitions. 

3.6 Forward-backward asymmetries 

The forward-backward asymmetry in B — * V£ + £~ decays is defined by 

dT/dq 2 \Jo dq 2 d cos 9 J-i dq 2 d cos 9 J 

oc Re (C%\{q 2 ) - R ut e ±ia C^|(g 2 )) , (37) 

where the second line follows from (jllj) . Note that we have not performed an average 
over CP-conjugate decays for the forward-backward asymmetry. The exponential reads 
e -ia f Qr g decay, and e la for B decay. 
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Figure 7: The forward-backward asymmetry in a) B + — > p + £ + £~, b) B~ — > p~£ + £~, 
and c) the CP-averaged S — > p°£ + £~ decay. The solid (dashed) line shows the next-to- 
leading (leading) order result. The band represents the theoretical error due to hadronic 
uncertainties. 



The next-to-leading order prediction of the forward-backward asymmetry for the 
B — > K*£ + £~ decay has been discussed in detail in our previous paper [3]. For the 
b — > s transitions the term Cg^]_(q 2 ) is negligible, because the corresponding R ut is very 
small. Hence there is no difference between B and B decay, and the asymmetry zero is 
determined by the zero of the real part of Cg ±(q 2 )- In [3J we found that the next-to- 
leading order correction shifts the zero by 30%, but once this correction is included, a 
precise measurement of the location of the zero translates into a determination of the 
Wilson coefficient Cg with an accuracy of about 10%. Our updated result for the position 
of the forward-backward asymmetry zero reads 

q 2 [K*°] = 4.36±°f GeV 2 , q 2 [K* + ] = 4.15±°J? GeV 2 . (38) 

The small difference compared to [S] is due to the different treatment of form factors 
and the inclusion of isospin breaking power corrections in the present analysis. 

In case of B — > p£ + £~ decays there exists an important new contribution from 
Cg U ]_(q 2 ). As a consequence, the decays of B or £>, neutral or charged B mesons to 
p£ + £~ may show significantly different forward-backward asymmetries. When a is near 
90° as expected in the Standard Model, we may approximate e ia ~ % sin a, and therefore 
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the additional contribution to the forward-backward asymmetry involves approximately 



Rut sin a Im (c^ (q 2 



(39) 



i.e. the absorptive part of the amplitude from the up-sector. We see from Table El 
that large contributions to the absorptive part arise from the virtual corrections to the 
b — > d£ + £~ transition, Y^ u '(q 2 ) and and that, in particular, the two-loop correction 
C (1) calculated in E| is as large as the one-loop term Y ( - U \q 2 ). The other large 
contribution comes from weak annihilation and exists only for charged B decays. In 



Figure U\ we show the expected forward-backward asymmetries for B 



± 



P 3 



and 



the asymmetry for the CP- averaged B 
asymmetry zeros are determined as 



P 



decay rate. 10 The locations of the 



B { 



p°£ + t 



<7o 2 [GeV 2 



A oa +0.02 1 +0.45 I 

-0.02lcKM-0.39lhad 



B+ -> p+ £+£- 

o c fi +0.09 1 +0.48 I 
0.00 -0.12lcKM-0.57lhad 



B° -> p° £ + £~ 

a i i +0.02 1 +0.26 1 

-0.02lCKM-0.22lhad 

/i /I o +0.07 1 +0.36 1 
^■^°-0.06lcKM-0.34lhad 



(40) 



As in the case of B — > K*£ + £~, the location of the asymmetry zero is a measure of C$ 
(now in the 6 — > ti sector), but in addition information about the phase can be obtained 
through the interference with the tree-dominated amplitude in the up-sector. 



4 Conclusions 

Our analysis provides Standard Model expectations for the exclusive, radiative and elec- 
troweak penguin decays B — > and B — > V£ + £~ (V = K*,p,u), extending our 
previous work [3j to b — > c?7 and 6 — > d£ + £~ transitions. The theoretical framework is 
complete at next-to-leading order in the strong coupling (except for weak annihilation 
which is included only in leading order), and leading power in the heavy-quark expan- 
sion. We have also included l/m^ power corrections, mainly from weak annihilation, 
which are important for the isospin asymmetries. 

The observables related to B — ► /ry decays can provide interesting constraints on the 
CKM triangle as was also recently discussed in [I5|ITEj. The accuracy of these constraints 
is limited by hadronic uncertainties, mainly from the hadronic form factors describing 
the B — > p transition (for branching fractions), 5-^7 transitions/ weak annihilation (for 
the isospin asymmetry), and from the scale and charm quark mass uncertainty (for direct 
CP asymmetries) . Nevertheless, the simultaneous measurement of branching fractions, 
isospin- and CP asymmetries can be used for an independent determination of the apex 
(p, fj) of the CKM triangle, which is complementary to the standard fit. In particular, if 
B° — ► p°7 is not observed with a branching fraction near the current upper limit 4- 10~ 6 , 
a tension with the standard fit arises. 

10 The case of to is not significantly different from p°. 
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We compared our calculations for B — > K*£ + £~ with the first experimental results 
on these decays. The central value of the Belle data on the partially integrated lepton 
invariant mass spectrum is about a factor of 2 larger than the theoretical prediction. 
Although the discrepancy is not significant within current uncertainties, it will be inter- 
esting to see how it develops. If the problem is theoretical, this may result in the curious 
situation that the comparison with data seems to favour a larger form factor A , but a 
smaller tensor form factor 7\. The present experimental accuracy for the measurement 
of the forward-backward lepton asymmetry in this decay mode is not yet competitive 
with the theoretical one. In the future, however, we expect significant constraints on the 
Wilson coefficient Cg in b — > s£ + £~ transitions as already discussed in [3]. 

We performed the first next-to-leading order calculation of exclusive b — > d£ + £~ 
transitions. In the long run when detailed experimental information on these transitions 
will become available, the CKM matrix can be assumed to be known. Hence the various 
observables related to B — > p£ + £~ decays provide information on the modulus and 
phase of the Wilson coefficient Cg in b — > d£ + £~ transitions or related coefficients in 
extended flavour models, which is of interest in the context of new-physics scenarios with 
non-minimal flavour violation. The theoretical uncertainty related to hadronic input is 
currently sizeable, but we may look forward to improved determinations of B —>■ p form 
factors, and perhaps, the mechanism of weak annihilation, on the time scales of interest. 
These improvements will add B — > p£ + £r to the list of rare processes that play an 
essential role in the program to uncover the origin of flavour violation. 
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A The amplitudes 

In this appendix we present the generalization of the amplitudes T^ 1 (a =_L, ||) given 
in [3| for B — > K*£ + £~ to the case of p or uj in the final state, and the new amplitudes 
as defined in Q. In 3| we used a set of Wilson coefficients denoted by C{. The 
exact relation between the Ci and the Cj given in Table ^ can be found in the appendix 
of 0- 

The coefficient functions CW and Tq ± defined in (jlj have the expansions 

= + + (41) 
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T® ± (u,u) = lfi>(ti,u;) + ^T^{u,u) + .... (42) 

In the following we give the expressions for the coefficients and T®±. The strong cou- 
pling is evaluated at the scale [i ~ nib in (jH]) and at /i/ = (0.5 GeV/i) 1 / 2 ~ (mbAQco) 1 ^ 2 
in (|42|) . which corresponds to the typical virtualities in the two terms. In contrast to our 
earlier analysis, we always evaluate the Wilson coefficients at \i nib, since the running 
of the four-quark operators ends at this scale. 



A.l Form- fact or term 

The coefficients Cj^'^ follow from Eqs. (12,14) in 3 . The corresponding expressions for 
(j(p,u) are obtained by the replacements C| ff — > and 

Y( s ) - y («) ( s ) = (to x + C 2 ) [h{s, m c ) -h(s,0)}, (43) 



where h(s,m q ) is given in Eq. (11) of 

The first-order corrections C^ 1 '^ are divided into a "factorizable" and a "non-factor- 
izable" term according to C^' 1 ' = Cj/' 1 ^ + C^ n ^' l \ The factorizable correction reads 



C (M = C eS I ln _ L + AM \ ; (44) 



Cf ') = - C f (in ^ + 2L + Am| (45) 

with L defined in Eq. (36) of [3 . AM depends on the mass renormalization convention for 
the overall factor nib in (J3J), such that AM = in the MS scheme, AM = 3 ln(m 2 /yU 2 ) — 
4(1 — u,f /nib) hi the PS scheme (our choice) and AM = 31n(mf/yU 2 ) — 4 in the pole 
mass scheme. Note that the expression for Cm differs from Eq. (35) of joj due to the 
different convention for the longitudinal "soft" form factor (jSJ). Furthermore, when 7]y 
is defined with Ti rather than the corresponding C^/'^' is given by ()44|) with the term 
— L omitted. The factorizable corrections from H^g follow again from Cj S — > 0, which 
leads to C ( /' u) = C\(' u) = 0. 

The non-factorizable corrections Cjf 1 ^^ are given in Eqs. (37,38) of [3], and make use 
of the result from The corresponding expressions for C^ 1 ^'^ are obtained by the 
replacements F 8 (7 < 9) - and -> Fg 9) + with given in 0. 



A. 2 Spectator scattering 

The longitudinal amplitude receives a leading-order contribution from a weak annihila- 
tion topology, where the photon couples to the spectator quark in the B meson [3] . It is 
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given by (e q denotes the charge of the spectator quark) 



rr(o,t)f \ M b uj AM b 34 

7ii (u, uj) = — e — ; C„ , 46 

where 

Cf = C 3 + ^(C 4 + 12C 5 + 16C 6 ) + 6 [1 + (-1)'] 5 gd (C 3 + 10C 5 ) (47) 

with J = 0, 1 for uj and p mesons, respectively. The corresponding expression for Jj|°^ 
is obtained by the replacement C^ 4 — > — C 12 with 

Cf = 3 V C 2 + 5 gd (J Ct + C 2 ) . (48) 

The terms proportional to 5 ? <j in the previous two equations do not appear in B — ► 
K*£ + £~ transitions. (Analogous terms with 5 q d — > <5 gs would appear in _B S — > <p£ + £~ .) 
There is no leading-order contribution from spectator scattering to Tjj ^ and the trans- 
verse amplitudes. 

The first-order corrections are again divided into a "factorizable" and a "non-factor- 

izable" term, T^± = T^± + T^± ■ There is also a first-order correction to the annihi- 

( f t) (ft) 

lation mechanism discussed above, which is not yet known. T±± and Tj| _ are given in 
Eqs. (20,22) of 3J, while now 

T (/,t) _ r ea 4M b (dq) 

because of the different convention pertaining to the longitudinal "soft" form factor. The 
corresponding terms from all vanish, T^±^ = 0. Furthermore, when is defined 
with T\ rather than the corresponding T±Y is zero. 

The non-fact or izable amplitudes T*£±^ are given in Eqs. (23-26) of The new 
amplitudes read 



M B ( 1 \ 
(u, uj) = e u ^— (C 2 - -Cxj [t±(u, m c ) - tj_(u, 0)) 

(u, uj) = e u (c 2 - ~Ci) (t\\ (u, m c ) - t\\ (u, 0)) , 



rp(nf,u) 

irrib 

n{nf,u) 



Tii (u,uj) = e, 



(nf,u), . M BUJ QM b 



9 Mbuj — q 2 — ie 

x (c 2 - -C^ (h(uM 2 B + uq 2 , m c ) - h(uM 2 B + uq 2 , 0)) , (50) 



6 

i( 1^ u \u,uj) = 0, where the functions tj_\\(u,m) are defined in Eqs. (27,28) of j^j. 
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A. 3 Power-suppressed amplitudes 

Some 1 / m^-suppressed weak annihilation contributions play an important role in decays 
to charged p mesons, because they are enhanced by the large Wilson coefficient Ci. In 
addition, power corrections may provide the dominant source of isospin breaking, since, 
as can be seen from the above formulae, the transverse amplitude is independent of 
the charge of the spectator quark in the leading order of the heavy quark expansion. 
On the contrary, all these effects are present in the longitudinal amplitude already at 
leading power. We therefore neglect power corrections to the longitudinal amplitude 
and summarize here the relevant expressions for isospin-breaking power corrections to 
the transverse amplitude. 

Weak annihilation. Denoting the power-suppressed contributions to T± defined in 
by ATj_ l \ we find for the annihilation terms at order a° s (s = q 2 /M B ) 



AT, 



(*) 



" e <?^ 7T~ Cs+o (C4 + 3C5 + 4t/ 6 ) / du 



AT, 



(u) 



3 m b M B V 3 J Jo u + us 

2n 2 f B f\\ m v 34 

+ q 3 m b M B (1 - s) X B , + (q 2 ) q ' 

2vr 2 / B /|| m v 



-(12 

3 m b M B (1 - s) \ B , + (V 2 ) ^ q ' 



„ ,7 ,^ Cf, (51) 



which generalizes the corresponding results in (SUES]- The inverse moments of the B 
meson distribution amplitudes, X B ,±(q 2 ), are defined in Eq. (49,50) of [3]. 

Hard spectator scattering. The power-suppressed hard scattering terms at order a s 
read [HI ES] 



AT, (<) 



= e q ^ T^rr-lnCf ^ hX ± (q 2 /M 2 B ) 
hsa q 4vr N c m b M B \ 8 M B j vy 1 Bj 



+ 8fJ 1 du^^F v t \uM B + u q 2 ) 
Jo u + us 

['du [ U dv^-F^\uM 2 B + uq 2 )\. (52) 
Jo Jo V 



(1 - s) X B , + (q 2 ) 

The quark-loop function Fy\s) is denoted Fy(s) in [3*H] . where also the integral X±(s) 
can be found. For q 2 = this integral suffers from a logarithmic endpoint singularity 
as u — > 1. We treat this singularity, which signals a breakdown of factorization for the 
power corrections, with the same ad hoc cutoff as in [351- The corresponding expression 



for AT, H 



is obtained by the replacements Cg ff — > and Fy\s) — > Fy \s), where 



Ff'(s) = - [C 2 - -d ) [/i(s,m c ) - /i(a,0)]. (53) 
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The numerically largest power correction is ATl , because it comes with a large 

ann 

combination of Wilson coefficients C 12 ~ 3 (when q — u, i.e. for B^ 1 decay). In this weak 
annihilation effect the photon is emitted from the spectator quark. Therefore the matrix 
element of factorizes at leading order in a s into 

(^(q,fi)p-(p',e*)\0^\B-(p)) = -if p m p e:(j*(q,^\ur(l- l5 )b\B-(p)) (54) 

assuming the decay is B~ — > p~£ + £~, for which the effect is most important. The 
B — ► 7* transition matrix element is dominated by hard scattering in the heavy quark 
limit, but power corrections (to weak annihilation, which is itself a power correction 
to the process of interest) may be significant, especially for q 2 ~ fIE\ I29j. where the 
hadronic structure of the photon may be resolved. To keep the discussion general, we 
parameterize the matrix element by 



flUl 



<7*(g,/iW(l - ls)b\B-(p)) = ig eui e q -^f^ V) [u^ %p n_ ff + 2<rf 

— 2£|| B7 \q 2 ) + corrections, (55) 

where the vectors n± are given by p' = En_ (E the energy of the p meson) and p = 
Mb(ji_+h + )/2. The form factors defined by this parameterization describe the transition 
of a B meson into a virtual photon of mass y/q 2 at zero momentum transfer p' 2 = 
(p — q) 2 = 0. In the heavy quark limit we find 

^ {q) -2X B+ (q 2 Y e " [q) -2X B ^q 2 Y (56) 

The "corrections" stand for additional Lorentz structures that vanish in the heavy quark 
limit, and also when the leading corrections from the hadronic structure of the photon 
are included. We may therefore adopt a more conservative treatment of the leading weak 
annihilation effects by expressing (|46|) (and, accordingly, T^^) in terms of £|| B7 \q 2 ) and 

(J51|) in terms of ^f 1 \q 2 )- The transverse form factor at q 2 = has been estimated with 
QCD sum rules |2*H1 129j. leading to a result not too different from the heavy quark limit 
(|56|). In our numerical analysis we use (|56|) to compute the central value of these form 
factors, but we assign a 50% theoretical error to this estimate. 
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